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Giant terahertz near-field enhancement by two-dimensional plasmons 
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We consider a periodically gated two-dimensional electron system, with a central gate finger biased indepen- 
dently creating a tunable plasmonic cavity in a planar plasmonic crystal. We demonstrate that the plasmons 
resonantly excited in the periodic plasmonic lattice by incident terahertz radiation can strongly pump the 
cavity plasmon modes leading to a deep sub- wavelength concentration of terahertz energy with giant electric 
field enhancement factor up to 10^. 
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Many material properties are sensitive to terahertz 
(THz) radiation, which makes THz spectroscopy an es- 
sential tool for the analysis of material microstructure, 
chemical composition and for the medical treatmentpii^. 
However, the diffraction limit of conventional optical sys- 
tems, i.e. their inability to focus light in deep subwave- 
length regime, challenges the study of individual micro- 
and nanostructures. 

Recent works have extensively discussed concentra- 
tion of the THz energy beyond the diffraction limit with 
metallic structures, including, nanoapertures^, tapers^^ 
and metal edge^i. In particular, in Refsj^i^ it was shown 
that tapered metallic rods and wedges can focus THz 
and mid-infrared radiation into a nanometer-size spots. 
In Ref^, THz energy concentration in a metallic nanoslit 
beyond the skin-depth limit was demonstrated with field 
enhancement factors up- to 10^ at 0.1 THz, neverthe- 
less, the effect was steeply decreasing with frequency in- 
crease (with factors below 100 at 1 THz). One of the 
key approaches for the design and study of the metal- 
lic nanostructures for THz photonics is making paral- 
lels with optical plasmonics. However, at THz frequen- 
cies the electromagnetic field is weekly localized near the 
metal surface and can be confined only between closely 
placed metallic surfaced or in dispersion-engineered ge- 
ometries^' Efficient coupling of THz radiation into the 
metallic nanocomponents demands comprehensive an- 
tenna design, which, together with the lack in tunabil- 
ity, makes the development of functional integrated THz 
photonic devices based on metallic nanostructures a chal- 
lenging task. 

Employing structures supporting intrinsic two- 
dimensional (2D) plasmons^ i.e. waves of free electron 
density physically localized in a plane, at THz frequen- 
cies, may resolve these problems and change significantly 
the design principles of THz photonic devices. Among 
such structures are semiconductor quantum wells. 



supporting 2D electron gases, and recently discovered 
graphene^^, a natural 2D material. The potential 
prospectives of graphene for THz photonics are cur- 
rently an active research fieldi^r— . Contrary to metal 
structures, the properties of such materials can be 
easily tuned by the applied electric potentialii. The 
latter effect was utilized for the development of compact 
tunable THz sources and detectors based on the field- 
effect transistor (FET) structures'^. In RefsJ^iii the 
enhanced plasmonic response in a periodic grating gate 
FET was demonstrated, and in RefJ^ a substantional 
increase in responsivity of such detector combined with 
a micro-bolometer element created by a potential barrier 
via splitting grating-gate was achieved. However, to the 
best of our knowledge, the THz field concentration and 
its enhancement in such 2D plasmonic structures has 
not been discussed yet. 

THz 
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FIG. 1. (Color online) Schematic of the split-grating-gate 
field effect transistor, where the central gate contact is biased 
individually. 



In this letter, we propose a novel mechanism for tunable 
THz near-field enhancement at a defect in a planar plas- 
monic crystal. We study plasmon excitation in the FET 
structure with a split-grating gate, where the central gate 
finger is biased independently from the rest of the grating 
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gate. We show that such configuration creates a tunable 
plasmonic cavity in the a planar plasmonic crystal lat- 
tice. We analyze the excitation of the lattice and the 
cavity plasmons, and demonstrate that their interaction 
leads to a dramatic THz near-field enhancement. 

We assume that a THz plane wave, with electric field 
perpendicular to the grating gate fingers, E 1 1 x, is in- 
cident normally upon the structure as shown schemati- 
cally in a Fig. [TJ A periodic metallic grating with period 
L = 4/im (2/im metal, 2/im gap) is placed on the top of 
GaAs substrate with dielectric permittivity e = 12.8. We 
suppose that 200nm below the wafer surface a 2D elec- 
tron gas with the electron density Nq = 2.5 x lO^^cm"^ 
is formed. 




FIG. 2. (Color online) (a) Plasmonic cavity resonance fre- 
quency as a function of 2D electron density under the defect 
gate for four cavity modes. Dashed line shows the position of 
the plasmonic lattice resonance. Dots mark the positions of 
the 3-rd and 4-th order cavity resonances shown in panel (b) . 

(b) Cavity excitation factor spectrum for the cavity depletion 
ratio Nc/No = 0.12. Numbers indicate the resonance orders. 

(c) ,(d) The mode profiles at the lattice resonance (cj = 0.336 
THz) and at the 4th order cavity resonance (cj = 0.388 THz), 
respectively. 



Negative electric potential applied to the central gate 
finger depletes the 2D electron channel under this partic- 
ular finger. Therefore, the biased area acts as a tunable 
plasmonic cavity (defect) in a periodic plasmonic lattice. 
Note that, the grating gate, if biased, also creates pe- 
riodic modulation of electron concentration in the 2D 
electron channel, and hence the frequency of the plasmon 
resonance of the plasmonic lattice can be also changed. 
Here, without loss of generality, we consider that the 
grating gate is not biased and the 2D electron density 
in the grating gated regions of the channel, A^o, remains 
homogeneous. 

When the frequency of the incident THz wave, cj, is 
far away from the lattice resonances, the lattice plas- 
mons can not be excited and the 2D channel response is 
governed mainly by the Drude background conductivity. 
In this case, neglecting the excitation of the lattice plas- 
mons, the plasmonic cavity can be fairly well modeled by 
a homogeneous 2D electron system gated only by a single 
gate finger. This system was analyzed in details in Ref^i^, 
where it was shown that the cavity resonance frequencies 
roughly satisfy a simple electrostatic formula: 



(/e) _ / Ncc'^h irk 

where k is the number of the resonance order, Nc is 2D 
electron density in the channel under the cavity gate, h 
is the gate-to-channel separation, and a is the width of 
the gate finger. 
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FIG. 3. (Color online) (a) Cavity excitation factor spectrum 
for different depletion ratios, Nc/Nq. Dashed line shows the 
position of the plasmonic lattice resonance, (b) Dispersions 
of the 3-rd and 4-th cavity resonances as functions of the 
the cavity depletion ratio. Solid curves show the dispersion 
of corresponding unperturbed cavity modes. Dashed straight 
line marks the position of the plasmonic lattice resonance, (c) 
The power stored in the unit cell of the plasmonic lattice with 
increasing the lateral distance from the cavity, (enhanced 
online) 

In Fig. [2fa), we plot the dispersion of the cavity reso- 
nance frequency defined by Eq.([T]) for four cavity modes 
with the variation of electron density in the biased area of 
the channel. In the same figure, we mark the position of 
the fundamental resonance of the plasmonic lattice. Ac- 
cording to Fig. [2fa) the dispersions of the higher-order 
cavity resonances cross the plasmonic lattice resonance 
dispersion with depletion of the cavity channel. In this 
case, we expect strong interaction between the cavity and 
the lattice plasmon modes. 

To prove our predictions, we perform comprehensive 
numerical simulations using the finite element method. 
We perform 2-D simulations, considering no variation of 
the field in the z— direction and searching for the field 
profiles in the (x-y) plane. We approximate the grat- 
ing gate by the sequence of infinitely thin perfect electric 
conductors. To model the structure, we use the super- 
cell approximation Ref^^. In this method, we embed our 
structure of interest into a large supercell periodically re- 
peated in the x— direction. We consider that the supercell 
contains the defect and 36 periods of the grating gate on 
either side of it. The electromagnetic response of 2D elec- 
tron gas in a local approximation can be described by the 
2D Drude conductivity, (72d{(^) = A/'(x)e^r/m*(l — jcjr), 
where N{x) is the 2D electron density {N{x) = Nc in the 
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FIG. 4. (Color online) (a) Normalized electric-field profile in 
the 2D electron channel and (b) maps of the electric-field 
distribution in the (x-y) plane for r = 10~^^s; (c),(d) the 
same as in (a) and (b) for r ■ 



10-^°s 



cavity and N{x) = A^o in the rest of the channel), r is 
the electron scattering time. 

We introduce the cavity excitation factor as a ra- 
tio between the THz power absorbed in the cavity 
and the energy flux in the incident THz wave, F = 
^Re{j^^^^^^(j\Ex\^dx)/ Pin- First, we demonstrate the 
vahdity of our approach for the prediction of the cav- 
ity resonance positions. Using Eq. ([1]) and Fig. [2fa), we 
choose the the depletion ratio Nc/Nq = 0.12 in such a 
way that the 3-rd and the 4-th cavity resonances are rel- 
atively close to the plasmonic lattice resonance. We plot 
spectrum of the cavity excitation factor for this case in 
Fig.[2fb). To identify the excited plasmon modes at each 
resonance, we analyze corresponding electric field pro- 
files, see Fig.[2jc,d). The analytically predicted positions 
of the lattice and the cavity resonances reasonably well 
coincide with the results of numerical simulations, see 
Fig. [2fa). The resonances of the cavity modes with the 
eigen-frequencies close to the plasmonic-lattice-resonance 
frequency become stronger due to the strong intermode 
interaction; whereas the cavity modes are negligibly weak 
far away from the lattice plasmon resonance. 

Moving further, we study the interaction between the 
lattice and the cavity plasmon modes. The cavity plas- 
mon resonances are shifted to lower frequencies for de- 
creasing depletion ratio, see Fig. [2l For depletion ra- 
tio Nc/Nq 0.1 the the 4-th order cavity resonance 
merges with the lattice plasmon resonance leading to the 
increase of the resonance amplitude, see Fig. [3fa) and 
corresponding media. For even smaller Nc/Nq ^ the reso- 
nances are separated again decreasing in amplitude, un- 
til the the 5-th cavity resonance approaches the lattice 
eigen- frequency. In Fig. [3jb), we plot the positions of 
the cavity resonances as functions of the depletion ra- 



tio showing the interaction between the cavity plasmon 
modes and the lattice plasmon mode. When the cavity 
and the lattice plasmon resonances approach each other, 
they experience strong interaction leading to the avoided 
crossing of corresponding resonances. 

In Fig.m^a) and (b), we plot the electric-field enhance- 
ment factor, |E/Eo| and the map of the electric-field dis- 
tribution in the (x-y) plane in the avoided crossing regime 
for Nc/No = 0.1 and = 0.332 THz. The electric field is 
greatly enhanced in the cavity up-to 150. The plasmon 
excitations in the lattice unit cells adjacent to the cav- 
ity are suppressed. With increasing the lateral distance 
from the cavity, the plasmon energy in the lattice cell sat- 
urates at the level corresponding to the lattice plasmon 
resonance, see Fig[3](c). This is due to the fact that the 
cavity mode is strongly pumped by the resonantly ex- 
cited lattice plasmons. The Q-factor of the resonances is 
about 20 30 for a room temperature electron scattering 
time r = 10~^^s. We have also studied THz near-field en- 
hancement for longer electron scattering time r = 10~^^s 
(electron mobility fi ^ 10^cm^/V)corresponding to cryo- 
genic temperatures. In this case, the Q-factor is much 
higher up to 3000 and hence the interaction between the 
cavity and the lattice plasmon modes is much stronger 
leading to a dramatic THz near-field enhancement factor 
|E/Eo| 0^ 10000, see Figs, g^c-d), which, to the best of 
our knowledge, is much greater than any of previously 
reported values^i^. Our calculations show that such dra- 
matic enhancement of THz near field by 2D plasmons 
can be achieved at higher THz frequencies (up to several 
THz) in the structure with a smaller grating-gate period 
and a sub-micron plasmonic cavity. 

In conclusion, we have studied THz field enhancement 
in the planar plasmonic crystal formed by a splt-grating- 
gate FET structure with a defect induced by an indi- 
vidually biased central finger. We have shown that the 
interaction between the cavity and the lattice plasmon 
modes can be electrically tuned by the voltage applied to 
the central gate finger. We predict unprecedented THz 
electric field enhancement factor up-to 10000 with a deep 
subwavelengh concentration of the THz electric field. We 
believe that such strong THz electric field enhancement 
can be utilized for THz nonlinear applications and for the 
design of the next generation of THz plasmonic devices. 
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